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Abstract. We present a detailed study of the cluster pair 
SL353 & SL349. This candidate binary cluster is located at 
the northwestern rim of the LMC bar Based on photomet- 
ric data we find that both clusters are coeval with an age of 
550±100 Myr We use the Call triplet in the spectra of indi- 
vidual red giants to derive radial velocities. Both components 
of the binary cluster candidate show very similar mean veloci- 
ties (« 274 ± 10 km s^i for SL 349 and w 279 ± 4 km s^^ for 
SL 353) while the field stars show lower velocities (« 240 ±19 
kms~^). These findings suggest a common origin of the two 
clusters from the same GMC. In this sense the cluster pair may 
constitute a true binary cluster. We furthermore investigate the 
stellar densities in and around the star clusters and compare 
them with isopleths created from artificial, interacting as well 
as non-interacting star clusters. Gravitational interaction leads 
to a distortion which can also be found in the observed pair. 

Key words: Magellanic Clouds - Hertzsprung-Russel (HR) 
and C-M diagrams - star clusters: SL 353 and SL 349 - Stars: 
kinematics 



1. Introduction 

The Magellanic Clouds offer the unique possibility to study star 
clusters in general and binary clusters in particular These two 
companion galaxies are close enough to resolve single stars, 
but distant enough to make the detection of close pairs of star 
clusters an easy task. About a decade ago, Bhatia & Hatzidim- 
itriou (119881), Hatzidimitriou & Bhatia (|1990|), and Bhatia et 



al. (1991) surveyed the Magellanic Clouds in order to cata- 



logue binary cluster candidates. To qualify as a binary cluster 
candidate, the maximum projected centre-to-centre separation 
of the components of a pair was chosen to be w lf3, which 
coiTesponds to « 19 pc in the LMC if a distance modulus of 
18.5 mag is adopted. A binary cluster with larger separation 
may become detached by the external tidal forces while shorter 
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separations may lead to mergers (Sugimoto & Makino 1989, 
Bhatia 1990| ). The number of chance-pairs of objects uniformly 
distributed in space can be estimated adopting a formula pre- 
sented by Page (1975): Roughly half of the pairs found may be 
explained by mere chance line-up. This suggests that at least a 
subset of them must be true binary clusters, i.e., clusters that are 
formed together and/or may interact or even be gravitationally 
bound. 

Depending on their masses and separations binary clusters 
will eventually merge or become detached. At one stage dur- 
ing the merger process the former binary cluster could have 
one single but elliptical core (Bhatia & McGillivray 1988). Re- 
cently deOliveira et al. (200C) performed numerical simula- 
tions of star cluster encounters which could represent a pos- 
sible scenario to explain the ellipticities found in several star 
clusters in the Magellanic Clouds. 

Star clusters form in giant molecular clouds (GMCs), but 
the details of cluster formation are not yet understood (Elme- 
green et al. 1999). Fujimoto & Kumai (1997) suggest that bi- 
nary or multiple star clusters form through strong, oblique col- 
lisions between massive gas clouds in high-velocity random 
motion, resulting in compressed sub-clouds revolving around 
each other. Since the components of a cluster pair formed to- 
gether, they should be coeval or at least have a small age dif- 
ference compatible with cluster formation time scales. Binary 
clusters are expected to form more easily in galaxies like the 
Magellanic Clouds, whereas in the Milky Way the required 
large-scale high-velocity random motions are lacking. Indeed 
only few binary open cluster candidates and no globular clus- 
ter pairs are known in our Galaxy (see for example Subrama- 
niam & Sagar 1999). However, binary globular clusters are not 
expected to survive the gravitational forces of the Milky Way 
(Surdin[l99l|). 



In case of tidal capture two clusters would be gravitation- 
ally bound, but age differences are likely. Encounters of clus- 
ters can be traced using isodensity maps (deOliveira et al. 
1998| , Leon et al. |1999| ). Though vandenBergh ( |1996| ) sug- 



gests that tidal capture becomes more probable in dwarf galax- 
ies like the Magellanic Clouds with s mall v elocity dispersion 
of the cluster system, Vallenari et al. (1998) estimated a clus- 
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ter encounter rate of dN/dt ^ 1 • (lO^yr) 
capture of young clusters very unlikely. 

The formation of low-mass star clusters tends to proceed 
hierarchically in large molecular complexes over several 10^ 
years (e.g., Efremov & Elmegreen 1998| ). Leon et al. (1999) 
suggest that in these groups the cluster encounter rate is higher 
and thus tidal capture is more likely: Binary clusters are not 
born together as a pair with similar ages but are formed later 
through gravitational capture. An observational test of this sce- 
nario would require the detection of evidence of tidal interac- 
tions between clusters, whose age differences need to be com- 
patible with the survival times of GMCs. 

Another binary cluster formation scenario is introduced by 
Theis (1998) and Ehlerova et al. (1997): Exploding supernovae 
close to the centre of a GMC sweep up the outer cloud mate- 
rial within a few Myrs and accumulate it in the shell. The large 
amount of matter makes the shell prone to gravitational frag- 
mentation and may eventually lead to the formation of many 
open cluster-like associations (Theis et al. 1997, Ehlerova et al. 
1997). In case of a dense ambient medium outside the cloud 



This makes tidal Table 1. Observing log 



or a very massive original molecular cloud the shell can be 
strongly decelerated resulting in a gravitationally bound sys- 
tem of fragments or stars which can recollapse. A galactic tidal 
field acting on this recollapsing shell can split it into two or 
more large clusters, thus forming coeval twin globulars (Theis 



1998 ). These clusters may stay together for a long time, though 
they are gravitationally unbound. The evolution of their spatial 
separation mainly depends on the shape of the shock front at 
the time of fragmentation. 

We are studying binary cluster candidates in the Magellanic 
Clouds to investigate whether the cluster pairs may be of com- 
mon origin and if they show evidence for interaction. While it 
is so far impossible to measure true, deprojected distances be- 
tween apparent binary clusters, an analysis of their properties 
can give clues to a possible common origin. 

The bina ry clu ster candidate SL353 (or BRHT33b, see 
Bhatia et al. |l99l| ) and SL 349 (BRHT 33a) is located in the 
outer west ern par t of the LMC bar. Based on integrated colours, 
Bica et al. ( |l996| ) suggest that SL 353 & SL 349 are coeval clus- 
ters of SWB type V (Searle et al. 1980 1 which is in agreement 
with the findings of Vallenari et al. ( |l998| ). Leon et al. ( |l999| ) 
suggest that the two clusters may even be physically connected. 
Based on ages for a large sample of star clusters, derived on the 
base of integrated colours, Bica et al. ( 1996 ) propose an age 
gradient of the LMC bar Younger clusters are predominantly 
found in the eastern part, while older clusters of SWB type III 
and higher are concentrated to its western end. 

This paper is organized as follows. In Sect. ^ we describe 
the photometric data in general. Stellar density maps are pre- 
sented in Sect. || The following section describes the colour 
magnitude diagrams (CMD) for the components of the cluster 
pair Ages for each cluster are derived and compared with pre- 
vious studies. The spectroscopic data are described in Sect. ^ 
and radial velocities are derived in Sect. ^. A comparison of 
the observed isopleths with simulated ones for interacting and 
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Table 2. Transformation coefficients 



Filter 


Ci 


Zi 








[mag] 


[mag] 


Gunn g 


-0.060 ± 0.006 


1.685 ± 0.030 


0.206 ±0.019 


Gunn i 


0.051 ± 0.003 


2.439 ±0.016 


0.081 ±0.010 



non-interacting star clusters is given in Sect. ^. In the last Sect. 
we summarize and discuss the results. 



2. Photometric observations and data reduction 

The data were obtained on March 22, 1994, with EFOSC 2 
at the ESO/MPI 2.2 m telescope at La Silla. A 1024 x 1024 
coated Thomson 31 156 chip (ESO #19) was used with a pixel 
scale of 0'.'34 resulting in a field of view of 5.' 8 x 5. '8. These data 
were obtained with the Gunn g (which is similar to Washington 
M) and Gunn i (corresponding to Washington T2) filters used 
at the 2.2 m telescope. An observing log is given in Table |l[ 
During standard image reduction with MIDAS we noticed flat- 
fielding problems near the edges of the exposures. Thus, each 
image was cut out resulting in a field of view of 4.'3 x 3f5. An 
image of the binary cluster candidate is shown in Fig. |l]. 

Profile fitting photometry was carried out with 
DAOPHOT II (Stetson |l99l| ) running under MIDAS. 

The photometry was transform ed us ing the standard fields 
around SA 1 10 and SA 98 (Geisler 199C ) observed in the same 
night. 

We applied the following transformation relations: 



g-G^ZG + aG-X + CG-{G-I) 
i - I = zi + aj ■ X + d ■ {G - I), 

where X is the mean airmass during observation, capital letters 
represent standard magnitudes and colours, and lower-case let- 
ters denote instrumental magnitudes after normalizing to an ex- 
posure time of 1 sec. The resulting colour terms Ci, zero points 
Zi and atmospheric extinction coefficients a; are listed in table 



3. Stellar density maps 

The cluster pair SL 353 & SL 349 (see Fig. |l]) is located in the 
inner northwestern part of the LMC at the outer rim of the bar 
The projected distance between the clusters' centres is l'.2A. 
Assu ming a distance modulus of m — M — 18.5 mag (Wester- 
lund 1997) t his co rresponds to 18.1 pc. According to Sugimoto 
& Makino (1989) this is close to the maximum separation at 
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Table 3. Number of foreground stars towards the LMC calcu- 



Fig. 1. Gunn i-image of the cluster pair SL353 & SL349. 
North is up and east to the left. The field of view of this image 
is 4.'3 X 3f5. These clusters are located at the outer northwest- 
ern rim of the LMC bar. For the stars marked with circles also 
spectroscopic data were obtained (see Sect. Bb 




Fig. 2. Star density map of SL353 & SL349. The star den- 
sity between the cluster pair is enhanced. Note the distortion of 
SL349 towards the bigger cluster SL353. Similar distortions 
are seen around the more massive cluster 



which a binary cluster is still stable and will not be detached by 
the outer tidal forces of the LMC. 

We investigate the stellar density distribution in and around 
the clusters by counting the number of stars in square cells of 
20 pixels length (corresponding to 6'.'8 or 1 .7 pc). To make den- 
sity structures and thus possible signs of interaction better vis- 
ible we applied a 3 x 3 average filter for image smoothing. 
The pr ocedure is described in more detail in Dieball & Grebel 
(|1998|). 



lated from the data of Ratnatunga & Bahcall (1985), scaled to 
our field of view of 4.'3 x 3f5 
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The stellar density map is presented in Fig. Q. The star den- 
sity between the two clusters is enhanced compared to the sur- 
rounding field density. The enhancement is about 3 a above the 
field background and thus significant at a « 99 %-level. In ad- 
dition, the smaller cluster SL 349 seems to be warped towards 
SL353. Also the density profile of SL353 shows distortions: 
we see smaller arcs of enhanced density reaching to the north 
and south, and also to the northeast. Whether these features 
might be tidal tails due to an interaction with the LMC tidal 
field or even due to an interaction w ith the smaller companion 
SL 349, as suggested by Leon et al. (1999), is hard to decide. 



4. Colour Magnitude Diagrams and Isochrone fitting 

To derive the CMDs for each cluster we cut out a circular area 
centred around the optical centre of each cluster. The radius 
was chosen to be 100 pixels (corresponding to 34" or 8.3 pc) 
for SL 349, and 120 pixels (corresponding to 40'.'8 or 9.9 pc) for 
SL353. An investigation of the star density plots suggests that 
most cluster stars are located inside the chosen areas. The clus- 
ter CMDs were compared with the CMD of the surrounding 
field and a statistical field star subtraction was applied to clear 
the cluster CMDs from contaminating field stars which might 
affect the age determination. To take into account the different 
area sizes of the regions covered by the clusters and the sur- 
rounding field we weighted the field star subtraction with field 
sizes. 

We derived ages for the clusters by comparing the CMDs 
with isochrones which are based on the stellar models of the 
Geneva group (Schaerer et al. 



1993). The isochrones were 



transformed to the Washington system by Roberts & Grebel 



(1994) 



A distance modulus of 18.5 mag (Westerlund 1997) was 
adopted. We assumed a metallicity of Z = 0.008 corresponding 
to [Fe /H] « —0.3 dex which was found by various authors for 
the young field star population (e.g., Russell & Bessell 1989, 



Jasniewicz 1992) 



Luck & Lambert 1992, Russell & Dopita 1992 Thevenin & 



Galactic field stars may contaminate our observed area. 
In order to assess whether foreground stars might affect our 
age determination we compare our CMDs with the number of 
galactic field stars towards the LMC estimated by Ratnatunga 
& Bahcall (1985). In Table ^ we present their counts scaled to 
our field of view of 4.'3 x 3.'5. As can be seen, the number of 
expected foreground stars is small, and we do not expect an 
influence on our age determination. 
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Fig. 3. Colour magnitude diagrams of the cluster SL 353 

(upper left) and its companion SL349 (upper right). The 
data points represent all stars which are located in an an- 
nulus of 120 pixels corresponding to 40'.'8 around SL 353, 
and of 100 pixels corresponding to 34" around SL349, 
respectively. The open circles represent stars that remain 
after statistical field star subtraction while the small dots 
represent field stars contaminating the cluster area. Small 
open circles denote stars lying inside a smaller annulus 
(17" for SL 353 and 13'.'6 for SL 349) where crowding is 
most severe. The crosses mark stars for which also radial 
velocities could be derived; small crosses denote stars 
with velocities that are concordant with field velocities. 
The colour of the stars redwards of the main sequence 
may be caused by the crowding in the clusters' centres. 
Several possible isochrones are overplotted on the CMDs. 
Due to crowding a definitive age determination is diffi- 
cult. We adopt an age of 550 ± 100 Myr for both clusters. 
Note the difference in age when fitting the main sequence 
turnoff or the red clump. The lower diagram shows the 
CMD of the surrounding field populations which com- 
prise a mixture of ages. Stars with T2 > 18 mag and 
red giants are represented as smaller dots to keep the 
isochrones recognizable. See text for more details 
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The CMDs of the two clusters and the surrounding field 
and the isochrone fitting are described in more detail in the 
following. 

SL353: The CMD of this cluster is shown in Fig. ^ upper 
left plot. Open and filled data points represent all stars which 
are located inside a radius of 40'.'8. After statistical field star 
subtraction only the open circles remain. The CMD shows a 
wide blue main sequence. The width of the main sequence is 
caused in part by photometric errors (average seeing 1'.'4) and 
crowding. Also a pronounced red clump of He-core burning 
stars and some red giants can be seen. Most red clump stars 
remain after field star subtraction, an indication of the inter- 
mediate age of this cluster The ratio of red clump stars in the 
cluster and in the surrounding field is approximately 3:1. The 
small open circles denote the stars that can be found inside an 
annulus of 17" where crowding is most severe. The width of 
the main sequence and the scatter of the red giants is caused 
in part by this crowding. Overplotted on the CMD are several 
possible isochrones. An isochrone resulting in an age of 500 
Myr (upper line) gives a good fit to the apparent main sequence 
turnoff at T2 « 18.5 mag. However, the evolved stars are not 
well represented by this isochrone. Fitting the He-core burning 
giants results in higher ages of 630 Myr (middle line) or 800 
Myr (lower line) but also underestimates the luminosity of the 
main sequence turnoff. Since the fit to the main sequence pro- 
vides a lower age limit we adopt a mean age of 550±100 Myr 
for SL 353. 

Such a discrepancy between the main sequence turnoff 
and red giants luminosities was noticed in the CMDs of var- 
ious other LMC star clusters, two prominent examples are 
NGC 1866 and NGC 1850 (see Br ocato et al. |1989| L attanzio 
et al. |1991[ , Vallenail et al. |1994b| , Brocato et al. |1994| )- Unre- 



solved binary stars which could increase the luminosity of the 
main sequence turnoff were proposed as one explanation. Other 
possibilities include blue stragglers and/or rotation (Grebel et 
al. 1996| ). Vallenari et al. (1994a) point out that no direct ev- 
idence for a significant population of binary stars could be 
found, though the detection of close binary stars in the Mag- 
ellanic Clouds is difficult if not impossible. Lattanzio et al. 
( 1991 ) was able to solve the problem by taking into account 
unresolved binaries as well as semiconvection. Their simulated 
CMDs then gave a good representation of the observed ones. 

SL349: Fig. ||, upper right plot, presents the CMD of 
SL 349. Again, large and small open circles denote stars that re- 
main after field star subtraction. SL 349 is the smaller one of the 
cluster pair, as can also be seen from the sparse main sequence 
and red giant clump. The small, open circles represent the stars 
located inside an inner, crowded annulus of 13'.' 6. Fitting an 
isochrone to the apparent main sequence turnoff (T2 k, 18.5 
mag) results in an age of 500 Myr (upper line). However, note 
again the pronounced red clump that is present in the CMD. 
We find approximately twice the amount of red clump stars in 
the cluster than in the field, scaled to the same area. Fitting 
isochrones to the core-He-burning red giants results in some- 
what higher ages than isochrone fits that take into account the 
apparent main sequence turnoff, namely 630 Myr (middle line) 



or 800 Myr (lower line). We adopt the same age of 550±100 
Myr for SL 349 that we already found for SL 353. 

Our findings are in agreement with former age determina- 
tions: Vallenari et al. (1998) found from CMDs a very simi- 
lar age of 500±100 Myr for both clusters. Based on integrated 
colours, Bica et al. (1996) suggest that the clusters are of some- 
what older SWB type V (800 - 2000 Myr). Studying the stellar 
content of proposed binary clusters in the LMC, Kontizas et 
al. (1989) found the older cluster pairs in their sample, among 



them SL353 & SL349, younger than 600 Myr, which agrees 



with our findings and those of Vallenari et al. (1998) 



The surrounding field: The CMD of the surrounding field 
is shown in the lower diagram of Fig. || A mixture of popula- 
tions of different ages can be found around the cluster pair. The 
CMD shows a bright blue main sequence and a few red super- 
giants which represent the young population. The intermediate 
age populations show up through the pronounced red clump of 
He-core burning stars and the red giants. Distinct young popu- 
lations cannot be distinguished, but the overplotted isochrones 
are supported by corresponding supergiants. Fainter main se- 
quence stars and red giants are plotted with smaller data points 
to keep the isochrones recognizable. The dotted isochrone fits 
to the brightest main sequence stars and to some supergiants, 
resulting in an age of 32 Myr. Also the 63 Myr isochrone (short 
dashed line) gives a good representation of the few red super- 
giants. We see again a few bright supergiants at T2 « 16.5 
mag which can be fitted with a 200 Myr isochrone (long dashed 
line). Between 200 Myr and 400 Myr (solid line) the star den- 
sity along the subgiant branch seems to be lower, indicating a 
possible decrease in the field star formation rate, but increases 
again along and below the 400 Myr isochrone (solid line). 

Since the main sequences in the CMDs show a large scat- 
ter, it is impossible to derive a reddening via isochrone fitting. 
Thus, we adopt a reddening of Eb-v — 0.1 mag, c orresp ond- 
ing to Em-T2 — 0.15 mag (see Grebel & Roberts 1995, their 
Table 5, for the transformation of extinctions in different filter 
systems), as suggested from the reddening maps of Schwering 
& Israel (|1991|). 



5. Spectroscopic data and their reduction 

In order to derive radial velocities, multi-object spectra (MOS) 
for 22 stars in and around the star clusters SL535 & SL349 
were obtained on February 09, 1995, using EMMI RILD at the 
ESO/NTT at La Silla. A TEK 2048 x 2048 chip (ESO #36) was 
used with a pixel scale of 0'.'27. The resulting field of view is 
9.'2 X 8.'6. Slit masks are prepared at the NTT during observa- 
tion, the punch field covers 5' x 8'. Only 19 stars out of 22 for 
which spectra were obtained are located inside the smaller field 
of view of the photometric data (4f3 x 3. '5). They are marked 
in Fig. [l|. The spectra were taken with grism #6 which covers a 
wavelength range of 6000 to 8300 A when the slits are centred 
on the CCD chip, and provides a dispersion of 1.2 A per pixel. 
The brightest stars in both clusters are red giants (see Sect. 
so we chose to obtain spectra of red giants at the Ca II triplet 
at A = 8498 A, 8662 A and 8542 A (see e.g. Olszewski et al. 
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Table 4. Observing log of the spectroscopic data Table 5. Radial velocities of all 22 programme stars 



Object 


Q2000 


^2000 


Exp.time 


velocity 




l^hm Sj 


r'"] 


[s] 


[kms-i] 


SL 353/349 


05 16 00 


-68 56 00.0 


3600 




HD 101266 


11 38 50.73 


-45 21 45.3 


30 


+20.6 


HD 111417 


12 49 31.82 


-45 49 32.8 


30 


-16.0 


HD 120223 


13 49 06.56 


-43 44 00.1 


30 


-24.1 



1991). The wavelength coverage of the MOS depends on the 
position of the slits with respect to the centre of the CCD. In 
order to reach the Call lines the X-position of the slits was 
shifted about « 2. '5 or « 530 pixels resulting in a wavelength 
coverage of 6500 to 9000 A. The slit width was « 1" (cor- 
responding to 4 pixels). In the same night the radial velocity 
standard stars HD 101266, HD 111417, and HD 120223 were 
observed using the same slit mask. An observing log is given 
in Tabl e |4|. T he velocities of the standard stars were taken from 
Evans dl967|). 



Data reduction was carried out using IRAF's twodspec 
and onedspec packages. We mainly followed the procedure 
for mu ltislit spectroscopic reductions described in Ellingson 
( 1989 ). Wavelength calibration was performed using He-Ar 



spectra obtained in the same night with the same slit mask. 



star no. 


velocity 


T2 


M - T2 


comment 




[kms^^] 


[mag] 


[mag] 




1 


313.7±3.1 





, 


variable/binary? 


2 


265.5±2.5 








field star 


3 


268.5±3.7 


14.52 


1.61 


field star 


4 


231.1±2.4 


16.62 


1.48 


field star 


5 


232.7±2.9 


16.31 


1.12 


field star 


6 


275.2±2.6 


16.67 


1.70 


cluster star 


7 


248.7±3.2 


16.57 


1.62 


field star 


8 


283.9±8.2 


16.97 


1.27 


cluster star 


9 


260.2±5.0 


15.24 


1.58 


cluster star 


10 


277.6±9.8 


16.47 


1.36 


cluster star 


11 


279.4±3.4 


16.56 


1.52 


cluster star 


12 


275.3±2.7 


15.74 


1.08 


cluster star 


13 


284.0±2.5 


16.35 


1.35 


cluster star 


14 


280.4±3.5 


16.35 


1.97 


cluster star 


15 


235.1±4.2 


16.11 


1.44 


field star 


16 


318.5±2.7 


16.53 


0.78 


variable/binary? 


17 


275.2±3.6 


16.82 


1.15 


cluster star 


18 


238.3±3.6 


16.67 


1.80 


field star 


19 


206.9±2.7 


16.56 


1.33 


field star 


20 


14.4±3.2 


16.05 


1.00 


foregr. star 


21 


229.0±3.3 


15.32 


2.02 


field star 


22 


84.3±3.4 






foregr. star 



6. Radial velocities 

We derived the radial velocities for the 22 stars by cross-corre- 
lating each spectrum with the spectra of the three standard stars 
listed in Table 0. An example is given in Fig.Q. The upper panel 
shows the spectrum of the Star no. 1 (which lies outside of the 
useful area of our imaging data) while the lower panel shows 
the spectrum of the radial velocity standard star HD 101266. 
The whole wavelength range covered by the grism is plotted. 
Strong atmospheric absorption lines are visible in both spectra. 
The sharp peak at A = 6867 A denotes the O2 B-band, and 
the even stronger double p eak at A — 7594 A belongs to the 
O2 A-band (Turnsheket al. 1985 ). Two H20-bands are visible, 
starting at A =7186 and 8164 A, respectively. 

Using the IRAF task fxcorr in the package rv we ap- 
plied a crosscorrelation between each programme and standard 
star spectrum. For this purpose we only need the Call lines. 
In Fig. Q the region used for the crosscorrelation is marked 
with dotted lines. Table || gives an overview of the radial ve- 
locities for all observed stars. The velocities listed are the er- 
ror weighted mean of the three velocities derived from the 
crosscoiTelation between the programme and the three standard 
stars. 

In Fig. H the spectra of all 22 programme stars are shown. 
We only plotted a wavelength range of 8350 A to 8750 A which 
includes the Call lines from which we derived the velocities. 
The two strongest lines of the Ca II triplet can be clearly seen in 
all spectra at w 8670 A and w 8550 A. The third and weakest 
line appears at « 8500 A. 

Star no. 8 is the faintest star of our sample and thus shows 
the noisiest spectrum (see Fig. ||), which is also responsible for 



the large eiTor of the velocity (see Table ||). Also star no. 10 
shows a somewhat noisy spectrum and has a higher cr. Stars 
no. 20 and 22 must be foreground stars since their velocities 
are small (14.4 km s~^ and 84.3 kms~^, respectively). The mi- 
nor displacement of the Ca II lines can be clearly seen in Fig. |[ 
All other stars have velocities which indicate that they are LMC 
members (approximatel y 200 - 320 kms~^, see e.g. Prevot et 
al. 1985 or Westerlund 1997 ). Stars no. 1 and 16 have veloc- 
ities higher than 300 kms^^. These two stars might be field 
stars having such a high velocity, or th ey may be variable or 
binary stars (see e.g. Fischer et al. 1993 ). 

The programme stars are located in the field as well as in 
both clusters and in between the two clusters. It is striking that 
most stars located inside the cluster areas have higher velocities 
while (almost) all field stars show lower velocities (see Table 
||). The components of our proposed binary cluster have very 
similar mean radial velocities of « 274±10kms^^ forSL349 
(stars no. 6, 8, 9, 10) and « 279 ± 4 kms^^ for SL353 (stars 
no. 11, 12, 13, 14, 17). The mean velocity for both clusters is 
« 277 ± 7 kms^^. In contrast the investigated field stars have 
velocities of « 240 ± 19 kms^^. The errors given for these 
mean velocities are the standard deviations. 

The expected internal velocity dispersion of a star cl uster i s 
of the order of only few kms^^ (see, e.g., Westerlund 1997 ). 
However, since the mean error of one single velocity measure- 
mentis w 6.6kms^^ we do not state about the internal velocity 
dispersion of SL 349 or SL 353. 
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Fig. 4. The upper panel shows the spectrum of the Star no. 1 while the lower panel shows the spectrum of the radial velocity 
standard star HD 101266. Strong night sky absorption lines are marked with arrows (only in the upper panel): The two sharp 
peaks at A = 6867 A and A — 7594 A denote to the O2 B-band and O2 A-band, the other two more extended bands starting at 
A =7186 and 8164 A, respectively, are caused by atmospheric water vapor. For the crosscorrelation only the wavelength range 
between the dotted lines was used. The Ca II triplet can be clearly seen in both spectra 



Table 6. Observational parameters used to create artificial star 
clusters 



cluster 


total mass 


number 


halfmass radius 


tidal radius 




[Mq] 


of stars 


[pc] 


[pc] 


SL353 


7300 


20100 


6.6 


23.0 


SL349 


4800 


13100 


3.3 


20.6 



7. N-body simulations 

Whether two clusters are interacting or not is hard to decide 
on the base of imaging data. According to deOliveira et al. 
( 1998 ) encounters of star clusters can be traced using isodensity 
maps. However, the field in which the clusters are embedded 
may likely influence the isopleths, and the distincti on be tween 
field and cluster stars is a difficult task (Leon et al. 1999). 



Assuming a Kroupa IMF (Kroupa et al. 1993) we estimate 
the total masses, the number of stars, and the half-mass and 
tidal radii of the star clusters SL 353 and SL 349. The observa- 
tional parameters can be found in Table ^. 

Based on a King model and the parameters listed in Table ^ 
two artificial star clusters were created which resemble the ob- 
served ones. Fig. ^ shows the density plot of the two artificial 
clusters placed at the same distance as the projected separa- 
tion between the cluster pair SL353 & SL349 (18.1 pc). The 



star density between the clusters, where the density profiles are 
overlapping, is enhanced. 

In order to get an idea whether an interacting cluster 
pair could be discerned from a mere overlap by chance in 
the line-of-sight, we performed N-body simulations using the 
GRAPE3a special purpose board in Kiel. The equations of mo- 
tion were integrated using the standard leapfrog scheme with 
a constant time step of 0.15 Myr. Additionally, we allowed for 
a Plummer softening with a softening length of 0.1 pc. Since 
we assume that the cluster pair SL 353 & SL 349 is a 'real' 
pair we chose a parabolic orbit as a limiting case to unbound 
pairs. In case of a 'true' gravitationally bound cluster pair, one 
would expect elliptic orbits and, thus, a series of repetitive en- 
counters which should result in even stronger signatures of an 
interaction or eventually already in a merger At the beginning 
we started with a well detached system at a separation of 100 
pc. The minimum distance (reached at t=0) was set to 12 pc. 
The simulation was stopped at separation of 18.1 pc, i.e. 13.4 
Myr after closest approch. For simplicity, the tidal field of the 
parent galaxy is not considered in this simulation. Fig. ^ shows 
a face-on view of both clusters. 

The isopleth of the interacting clusters, again at a distance 
of 18.1 pc, is presented in Fig. |^. Again, a region of enhanced 
star density connects both clusters. Furthermore the larger clus- 
ter shows a distortion while the smaller but more compact one 
remains mainly undisturbed in the isopleth. The density plot of 
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Fig. 5. Spectra of the stars no. 1 to 22. Only the wavelength range of 8350 A to 8750 A which includes the Ca II lines and which 
was used for the crosscorrelation is plotted. The Ca II triplet rest wavelenghts are at A = 8498 A, 8662 A and 8542 A. The Ca II 
lines for the stars no. 20 and 22 are less displaced compared to the other spectra which identifies them as galactic foreground 
stars. See text for more details 
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Fig. 6. N-body simulation of a close encounter of two star clusters which resemble our observed cluster pair The two clusters 
approach each other on elliptical orbits under the influence of their gravitational forces. The tidal field of the parent galaxy is not 
considered in this simulation 
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simulation is barely possible since the parameter space cannot 
be further restricted on the base of our data. 




Fig. 7. Density plot of two artificial star clusters without any 
interaction. The centre-to-centre separation between the clus- 
ters is the same as the projected distance between the observed 
pair SL353 & SL349 (18.1 pc). The star density between 
the clusters, where the density profiles are overlapping, is en- 
hanced 




Fig. 8. Isopleth of the two artificial clusters with interaction, 
the distance between the two clusters is the same as in Fig. 0. 
The region between the clusters shows enhanced star density. 
The larger cluster shows a distortion while the smaller but more 
compact one remains mainly undisturbed in the isopleth 



the observed cluster pair (Fig. ||) shows a similar, northeast to 
south distortion around the more massive cluster SL 353. How- 
ever, the observed northern distortion and the distortion of the 
smaller cluster SL 349 are not seen in Fig. ||. 

We want to emphasize that the parameters used to create 
the artificial clusters are based on our observations, but are only 
estimates and should not be taken as strict values. The simula- 
tions are meant for qualitative comparison and a more detailed 



8. Summary and conclusions 

Fitting isochrones based on the Geneva models (Schaerer et al. 



1993) to the CMDs, we find that the components of the double 



cluster SL 353 & SL349 are coeval within the accuracy of our 
data, and we derive an age of 550±100 Myr for both clusters. 

The clusters are sufficiently old that the Ca II triplet visible 
in the spectra of red giants could be used to derive radial veloc- 
ities. 22 stars in and around the two star clusters were investi- 
gated. Most stars located inside the cluster areas show similar 
velocities of « 277 ± 7 kms^^ whereas almost all field stars 
show lower velocities of « 240 ±19 kms^^. Two foreground 
stars were identified through their velocities which are too low 
to belong to the LMC population (14 kms^^ and 84 kms~^, 
respectively). 

Both components of the binary cluster candidate are of the 
same age and, furthermore, have very similar mean radial ve- 
locities of w 274±10kms"i forSL349 and « 279±4kms"i 
for SL 353. These findings support that both clusters may have 
formed at similar times and from the same GMC. In this sense 
they may constitute a true binary cluster. 

We investigate the stellar densities in and around the binary 
cluster SL 353 & SL 349 and find an enhanced density between 
the two clusters. The smaller cluster SL 349 shows a distortion 
towards SL353, and similar distortions can be seen around the 
more massive cluster. 



Vallenari et al. ( |1998| ) found a distortion of the isophotal 
contours, similar to our Fig. |[ and a twisting of the isopleths 
which they regard as a sign of interaction and physical connec- 
tion between the two clusters. It is remarkable that the age of 
this binary cluster is higher than the theoretical survival time of 
few 10^ yrs fo r physically connected cluster pairs suggested by 
Bhatia ( |1990| ). Leon et al. ( |1999| ) suggest that SL 349 & SL 353 
are part of a larger star cluster group in which interacting binary 
clusters are formed later via tidal capture, thus explaining an 
age larger than the theoretical survival time. However, another 
explanation could be that the survival time might be larger in 
the LMC bar where the tidal field might be weaker, as proposed 
by Elson et al. ( |1987| , their Fig. 13). 

Whether the distortion of the clusters is indeed a sign of 
possible interaction cannot be decided on the basis of our imag- 
ing data alone. 

Based on observational parameters we created two artificial 
star clusters. The isopleths of the non-interacting pair as well as 
of a simulated interacting pair of star clusters were compared 
with the observed density plot. The density plot of the artificial, 
interacting pair shows a distortion of the more massive cluster 
which can also be found in the observed isopleth, however, the 
observed distortion of the smaller cluster SL 349 is not seen in 
the artificial density plot. It seems likely that this cluster pair 
shows signs of interaction, however, this does not necessarily 
imply that both clusters are gravitationally bound. 
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